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We study coherent electronic transport in charge density tunable micro- devices patterned from 
thin films of the topological insulator (TI) Bi2Se3. The devices exhibit pronounced electric field 
effect, including ambipolar modulation of the resistance with an on/off ratio of 500%. We show that 
the weak antilocalization correction to conductance is sensitive to the number of coherently coupled 
channels, which in a TI includes the top and bottom surfaces and the bulk carriers. These are 
separated into coherently independent channels by the application of gate voltage and at elevated 
temperatures. Our results are consistent with a model where channel separation is determined by 
a competition between the phase coherence time and the surface-to-bulk scattering time. 
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Topological Insulators (TIs) are gapped bulk insula- 
tors with gapless Dirac surface states which have emerged 
as a new paradigm in the study of topological phases 
of matter Tl-based electronic devices are attractive 
as platforms for spintronic applications [H, and for de- 
tection of emergent properties such as Majorana excita- 
tions electron- hole condensates [1| and the topologi- 
cal magneto-electric effect (Hf . Most theoretical proposals 
envision an experimental geometry consisting of a planar 
TI device, where electrical current is carried by the sur- 
face states. 

Despite considerable recent evidence of TI surface 
states in ARPES [sl, and STM [8], their observation 
and manipulation in transport experiments remains dif- 
ficult: TI devices require a density-tunable surface state 
which is decoupled from the residual bulk carriers. Re- 
alistic devices, however, conduct through parallel chan- 
nels consisting of the top and bottom surfaces, and of 
bulk carriers which, due to unintentional doping, can 
account for a significant part of the conductance, limit 
the surface density tunability, and create an uncertainty 
in the surface-to-bulk coupling. It is therefore desirable 
to minimize the bulk contribution and to simultaneously 
investigate the various ways in which it is involved in 
electronic transport. Suppression of the bulk channel 
can be obtained by minimizing carrier density, as demon- 
strated in mm-size single-crystals [sl, and by fabrication 
of nanoscale devices, such as nanoribbons [10| and flakes 
exfoliated from single-crystals [lll-[l3| . The latter studies 
have also demonstrated that the surface carrier density 
can be tuned by electrostatic gating and have all detected 
an ambipolar modulation of conductance. 

An alternative route for fabrication of thin TI devices 
employs the growth of thin films in ultra-high vacuum 
chambers. Thin film growth offers fine control over ge- 
ometry and composition and a straightforward approach 
for growth of heterostructures. High quality Bi2Se3 thin 
films were studied by in-situ STM Q and ARPES 
Transport studies TgI, [l^ carried out on the same films 



^^electric 




20 (deg) 



FIG. 1 (color online): (a) Optical image of a gated Bi2Se3 
device, showing a Bi2Se3 Hall bar, evaporated contact pads, 
and top-gate. The gate electrode covers the entire area of 
the Hall bar. (b) Schematic of the device and measurement 
circuitry, showing the Bi2Se3 layer, cap layer, contact pads, 
dielectric, and top gate, (c) X-Ray Diffraction intensity of 
a Bi2Se3 thin film, showing the (0,0, 3n) family of diffraction 
lines indicating growth along the c axis. Some lines associated 
with the Si substrate are also visible. Inset: Image of Bi2Se3 
film grown on a Si wafer. 

are still dominated by bulk transport since the material 
is highly doped. Interestingly, most thin film studies re- 
port a pronounced weak antilocalization (WAL) feature, 
which is an indication of phase coherent transport. 

WAL and weak localization emerge from the correc- 
tion of coherent time-reversed closed paths to electronic 
transport [l^. They are sensitive to the competition be- 
tween the phase coherence time and other time-scales, 
and have been extensively employed as a probe for co- 
herent dynamics in solid state systems [l9|. In TIs both 
the bulk and surface states may contribute to coherent 
transport, and in both cases they should exhibit WAL: 
In the bulk, the strong spin-orbit (SO) coupling leads to 
random rotations of the spin orientation [20|, which on 
average result in destructive interference for backscatter- 
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FIG. 2 (color online): (a) Square resistance R vs. top-gate 
voltage Vtg measured on devices Al (continuous line), A2 
(dashed) and A3 (dash-dot) at T = 4K. Inset: Square con- 
ductance, C, for device Al. (b) Low field Hall coefficient 
Rh in 2D units. Rh is observed to cross over to positive 
values at Vtg ^ -7V. (c) Schematics of the band structure 
and spatial variation of bulk bands near the surface at dif- 
ferent gate voltages. In each subpanel, the left side shows 
the band energy vs. vertical position, and the right side is 
the energy-momentum dispersion at the surface. The bulk 
conduction band is marked by "Ec"; Valence band by "Ev"; 
and surface by "S". Right panel (hollow square): Vtg > 0. 
The bulk states near the surface are populated. Middle panel 
(solid square): Bulk states near the surface are at the deple- 
tion point (bottom of conduction band). Left panel (circle): 
Vtg < 0. Top surface state is gated to the Dirac point. The 
bulk bands bend to form a depletion layer near the surface, 
marked by vertical lines. This depletion region increases with 
increasing negative Vtg- The markers in panel (a) correspond 
to the band structure schemes shown in the bottom panels. 

ing paths, hence leading to enhanced conductance. On 
the chiral TI surface state, momentum is coupled to the 
spin degree of freedom, so time-reversed paths around a 
closed trajectory acquire a relative phase of tt [21]. 

Given the various parallel conduction channels, so far it 
was not clear why most of the recent TI transport studies 
report WAL corrections corresponding to a single coher- 
ent channel [IgI, [l^ , and it was suggested that one of the 
surfaces has a significantly higher coherence length than 
the other [l^- A few studies have found the number 
of channels to be tunable by a gate voltage [l3|, [iTI, |2^ , 



but the microscopic mechanism underlying this tunabil- 
ity was not thoroughly investigated. Here we study the 
WAL effect in charge density tunable Bi2Se3 devices. We 
find that the number of independent coherent channels is 
tunable by the applied electric field and by the tempera- 
ture. Our results indicate that bulk carriers play a crucial 
role in TI coherent transport, and that channel coupling 
is controlled by a competition between the phase coher- 
ence time and the surface-to-bulk coupling time. 

We fabricate low density Bi2Se3 devices by growing 



large area thin films on a Si(l 11) substrate and sub- 
sequent lithographic processing [22|. A device image 
and schematics are shown in Figs. 1(a) and (b). X-Ray 
diffraction (Fig. 1(c)) reveals sharp (0,0, 3n) peaks, con- 
firming that the films are c-axis oriented along the growth 
direction. The Bi2Se3 is typically capped in-situ by a 
sputtered layer of 4nm AI2O3. Additional 16nm of Hf02 
gate dielectric are deposited after device patterning. 

Figure 2 shows the electric field effect behavior of 
three representative devices (Devices Al,2,3), all pat- 
terned from the same 20nm thick capped Bi2Se3 wafer, 
and measured at T = 4K. Every resistance (conductance) 
shown is resistance (conductance)-per-square. For Vtg 
> 0, is modulated weakly by the gate, but begins rising 
sharply around Vtg = 0, and peaks at 5 - 6 kfl around 
Vtg = -5V. The inset shows the square conductance G of 
Device Al in units of e^/h. For this device G varies be- 
tween 5 and 25 e^//i , exhibiting a pronounced ambipolar 
modulation characteristic of Dirac dispersions, and indi- 
cating that the conductance is dominated by the top sur- 
face. This is supported by the Hall coefficient Rh^ which 
for a single channel is given by Rh = — 1/ne, where n is 
the density of charge carriers and e is the electron charge, 
but which is more complex for multi-channel systems |11|. 
Rh is strongly modulated by the gate and crosses over 
to positive values (Fig. 2(b)) at Vtg ^ -7V, indicating 
that the dominant Hall current-carrying population has 
changed from electrons to holes, as expected for Dirac 
systems with low doping. This is in stark contrast to 
previously reported studies on Bi2Se3 thin films, which 
are very strongly electron-doped, and do not exhibit such 
a strong carrier type modulation 16|, 17 1. Gating of the 
surface states shifts the surface energy bands vertically 
and hence is accompanied by bending of the bulk bands 
near the surface, as shown in Fig. 2(c). We associate 
the sharp change in the slope of R{Vtg) near Vtg = 
with the depletion of bulk carriers immediately near the 
surface. For Vtg < 0, the main effect of the gate is to 
change the TI surface state charge density, which results 
in a more rapid variation of conductance. However, any 
change of the surface charge density has to be accom- 
panied by a change in the width of the bulk depletion 
region. This depletion region is crucial for coherently 
separating the surface state from the bulk bands. 

We now turn to the electric field effect on coherent 
transport, studied by magnetoconductance. These data 
were taken on Device B, which is patterned from a dif- 
ferent wafer than the one reported above but has similar 
transport characteristics, with a well defined resistance 
maximum (Fig. 3(b)). Figure 3(a) shows R vs. perpen- 
dicular magnetic field, 5, at T = 0.3K and where WAL 
appears as a sharp suppression of resistance at low mag- 
netic field. We repeat the same measurement for different 
gate voltages, and find that the WAL feature also evolves 
with Vtg^ as shown in Fig. 3(c), where the change in con- 
ductance, AG{B) = G{B) - G(0), is plotted at Vtg = 
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FIG. 3 (color online): (a) magnetoconductance of Device B 
taken at Vtg — +6V showing a pronounced WAL suppression 
of resistance, (b) R{B = 0) vs. Vtg- R peaks at Vtg = -8V. 
(c) AG vs. B plotted for Vtg = +6V and -8V (marked by 
circles in (b)). The two plots are offset for clarity. The data 
are plotted in blue, and fits to WAL theory (Eq.[T]) are plotted 
in yellow, (d) Coherence length and prefactor a extracted 
from Eq. [T] for each Vtg- Both and a are modulated by 
the applied gate voltage. 



limit, TsB < charge carriers scatter between the bulk 
and surface states while maintaining phase coherence, ef- 
fectively becoming a single phase-coherent channel. We 
therefore interpret the gate-dependence of a as direct 
modulation of tsb via the formation of a depletion layer 
between the top surface and bulk carriers (see Fig. 2(c)), 
spatially separating them and suppressing the scattering 
probability between them. 

It is important to discuss the validity of Eq.[T]when the 
coupling between channels is tunable. In the fully decou- 
pled regime, each channel i has correction AGi which 
follows Eq. [1] (with alpha = -0.5) and depends on Z^^^, 
yielding for the total correction AGtot = AGi + AG2. 
Decomposing Eq. [T] into its logarithmic and digamma 
components [22|, we note that the latter approaches a 
constant value for B > Bq ^ 10 — 20mT. Since our data 
extend to a few T, where the change in AG is dominated 
by the logarithmic component, Eq. [T]is a very good ap- 
proximation to AGtot ^ with a = —1 and an effective 
coherence length l^^-^ = ^JT^li^^ as shown in [lH. As a 
consistency check, we can find the dataset where Eq. [T] 
yields a = — 1, and fit the data to AGtot [HI- We find 
the fit to agree very well with the data both for B > Bq 



+6V and Vtg = -8V. Both traces exhibit a correction ^ <^^-' ^^^^ coherence lengths /^i,2 = 135, 77nm. 

of AG ~ although the former is sharper in mag- 



netic field. The WAL correction depends on the phase- 
coherence time r0, or length Z^, which are related through 
1(f) = y^Dr^ {D being the diffusion constant). The data 
agree well with [23[: 
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where Bo = h/ (4e/^), i/j is the digamma function, and a is 
a prefactor which should be equal to -1/2 for a single co- 
herent channel. Eq.[T]is valid for massless Dirac fermions 
(i.e. the surface states) [2l[ and is also valid in the bulk 
where the spin-orbit scattering time rso is significantly 
shorter than [22| . The data (blue) agree well with the 
fit to Eq. [1] (yellow) over the entire measured range in 
both cases, indicating that WAL constitutes the entire 
correction to magnetoconductance. The fit contains two 
free parameters, the coherence length and the prefac- 
tor a, both of which change with the applied gate voltage 
(Fig. 3(d)). Here we focus on the gate-tunability of a, 
which changes from -0.7 to ~ -1, reflecting a change in 
the effective number of coherent channels. 

From the WAL point of view, the surface and bulk 
states can be regarded as independent phase coherent 
channels as long as the carriers in one channel lose co- 
herence before being scattered into the other. In this 
case, the majority of closed loops responsible for WAL 
will involve states from a single channel, and each chan- 
nel will exhibit its own WAL correction. This condition 
can be formulated as r^^ > t^, where tsb is the ef- 
fective surface to bulk scattering time. In the opposite 



When —0.5 > a > —1.0 the physical interpretation of 
the parameters extracted from Eq. [T] is not trivial and 
no theoretical model exists to predict the magnetocon- 
ductance in the crossover coupling region. The evident 
success of the fit indicates that the logarithmic correc- 
tion is robust, and a can be used phenomenologically as 
a measure for channel separation. 

We have so far encountered two effective regimes dif- 
fering in the degree of channel separation: (i) Vtg = -8V, 
where a ~ —1 = —(1/2 + 1/2), indicating a decoupling 
of the top surface from the rest of the system; (ii) Vtg = 
+6V, where a ~ —0.7, indicating that the top surface is 
only partially decoupled. In highly doped samples (De- 
vice C, Fig. 4(b)) we find a third regime, where a = -0.5, 
Rh = -3.7 n/T, and G = 57 eVh. This result is consis- 
tent with other WAL studies TgI, l^ and suggests that 
the surface and bulk channels in highly doped TIs are 
fully coupled. This is consistent with our model, since 
for Vtg > the surface and bulk states co-exist in space 
and I /tsb should depend on the momentum difference 
between the surface and bulk bands, which become closer 
as the density increases 0, [lH . 

The effective number of phase coherent channels re- 
flected in the magnitude of a depends on the ratio tsb 
to r^, which was controlled above by modulating tsb via 
the electric field effect. An independent control over this 
ratio can be obtained by varying the temperature, which 
directly controls through dephasing. Figure 4 shows 
the change in the WAL correction with temperature in 
the different mixing regimes discussed above. Fig. 4(a) 
shows AG{B) in regime (ii) from T = 1.5K to 47K, to- 
gether with the corresponding fits to Eq. [TJ Similar data 
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FIG. 4 (color online): (a) Magnetoconductance AG measured 
on low-density Device B at Vtg — OV for different temper- 
atures (from 1.5K to 47K). The curves are vertically shifted 
for clarity. Superimposed on each scan, in yellow, is the fit to 
Eq.[Tl (b) WAL traces and fits to Eq.[T] taken on low density 
Device B and high density Device C. Device B: = 175nm, 
a — -0.75; Device C: = 300nm, a — -0.5. (c) vs. T for 
the three regimes discussed in the text. Open circle: Device 
B at the resistance peak (i); Full circle: Device B at Vtg — 
OV (ii); Triangle: Device C (iii). (d) a vs. T, for all three 
regimes. 

and fits are obtained for regimes (i) and (iii). The corre- 
sponding temperature dependence of a and is shown 
in Figs. 4(c) and (d). a decreases with temperature for 
regimes (ii) (from -0.75 to ~ -1) and (i) (-1 to -1.15), 
suggesting increased channel separation with tempera- 
ture. This may seem surprising, since naively one expects 
larger channel mixing at higher temperatures. However, 
the behavior of a(T) is consistent with our model, and 
can be understood once the temperature dependence of 
1(f) is examined. In Fig. 4(c) we see that is strongly 
temperature dependent, and since Tfj^ ^ l'^^ this indicates 
that T(f) decreases rapidly as T is increased, which should 
result in a decrease in the t^Itsb ratio, and consequently 
an increase in the channel separation. The validity of this 
model requires that tsb changes slower than r^. This is 
expected in view of recent mobility data [11| suggest- 
ing that the impurity scattering rate, which should gov- 
ern TsB^ is indeed nearly temperature independent below 
40K. The results in regime (i), where ol becomes smaller 
than -1, suggest that, in addition to the top-bulk sepa- 
ration identified above, also the bottom surface could be 
decoupling from the bulk . At very high densities (Device 
C), this behavior is not observed, likely because tsb is 
much shorter than in the temperature range explored. 
The corresponding increase in a at high T could be con- 
sistent with coherent transport being dominated by the 
bulk, where WAL is suppressed when ^ tso- 



During the preparation of this manuscript we became 
aware of work by Chen et al., [2^ , which report a similar 
modulation of a, but associate it with changes in the 
coherence lengths of electron and hole channels. Such 
interpretation cannot explain our data, since our samples 
are thicker and the bottom surface density is not tunable, 
nor can it explain our observed temperature dependence 
nor the fact that a = -0.7 at positive gate voltages. 
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